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ABSTRACT 

The Sun and the undisturbed interstellar magnetic field and plasma velocity vectors 
{Bis,Vis) define a mirror symmetry plane of the flow at large heliospheric distances. 
We show that for the Bis direction defined by IBEX Ribbon center, the radial di¬ 
rection of Voyager 2 over the last decade, and the (thermal proton) plasma velocity 
measured by the spacecraft since 2010.5, are almost parallel to the {Bis,Vis)- plane, 
which coincides in practice with the Hydrogen Deflection Plane. These facts can be 
simply explained if approximate mirror symmetry is also maintained on the inner 
side of the heliopause. Such approximate symmetry is possible since the solar wind 
ram pressure is almost spherically symmetric and the plasma beta value in the inner 
heliosheath is high. In the proposed symmetry, the plasma flow speed measured by 
Voyager 2 in the inner heliosheath is expected to rotate more in the tran s verse than 
in the polar direction (explanation alternative to iMcComas fc Schwadronl (12014D ). in 
evident agreement with available spacecraft data (our Fig.l). 
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1 INTRODUCTION 

Both Voyager spacecraft continue their journey away from 
the Sun through the heliospheric interface region. Voyager 1 
( VI) has already crossed the heliopause (HP) dBurlaea et al.l 
120131 : [Stone et all 120131 : [Gurnett et al.ll2013h . the boundary 
that separates the heliospheric plasma from the interstellar 
one, and is now operating in the interstellar medium at a 
distance of 130 AU (early 2015) from the Sun. Voyager 2 
(V2), at 106.8 AU from the Sun, still penetrates the inner 
heliosheath region. 

At the time when Voyagers commenced their outward 
journey neither the strength nor direction of the interstellar 
magnetic field (ISMF, Bis) was known. In 2009, the Inter- 
stellar Boundary Expl orer ( IBEX) discovery of the ribbon 
llMcComas et al.l 120091 1 indicated the probable direction of 
the ISMF as the ribbon center (RC) llFunsten et al.ll2009l ). 


2 GEOMETRICAL COINCIDENCES 

It is noteworthy that each of the Voyagers trajectories 
happens to have a special orientation (instances of ge¬ 
ometrical coincidence) with respect to the direction de¬ 
fined by the RC. The first such case was discussed in 
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iGrvgorczuk. Czechowski fe Grzedzielskil (2014|). In that pa¬ 
per we pointed out that VI trajectory shares the same he¬ 
liographic latitude (^34.5°) with the RC. This led us to a 
simple explanation why the magnetic field directions mea¬ 
sured by VI on both sides of the HP were so similar. 

In the second case, to be discussed now, the V2 
spacecraft trajectory direction is close (within ~2° ) to 
the plane (the {Bis,Vis) plane) which contains the po¬ 
sition of the Sun and the directions of undisturbed in¬ 
terstellar plasma flow {Vis) and the IBEX RC ( Bis)- 
This plane was pointed out llGrygorczuk et ahl 1201 j) to 
be close to the hydrogen deflection plane (HDP, the 
plane defined by the interstell ar H a n d He inflow ve¬ 
locity vectors, fLallement et al.l (|2005l . 120101 )). In cur¬ 
rent heliospheric modeling, the HDP is generally iden¬ 
tified with the (Bis,Vis) plane (llzmod enov et ah 20051; 
Pogorelov fe Zank 2006; Ophcr et al. 2006: IPogorelov et al 


Heerikhuisen fe Pogorelovl 2011 


Katushkina et al 


If the supersonic solar wind were fully spherically sym¬ 
metric with no interplanetary magnetic field, the (Bis, Vis) 
plane would obviously had been the mirror symmetry plane 
of the problem. The symmetry is broken for a number of 
causes: fast/slow wind, angle and time dependent solar ac¬ 
tivity, solar rotation, Parker magnetic field, etc. However, 
if these effects are relatively weak, an approximate mirror 
symmetry should be discernible also in the inner heliosheath. 

In this letter we show that V2 measurements agree 
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All —sky Map in Heliographic Inertial Coordinates 



Figure 1 . Mollweide projection of the plasma velocity directions observed by V2 (years 2010 - 2014), as well as V2 locations. Shown are 
the yearly averages of raw data (black crosses) and six-month averages of the instrumentally corrected velocity directions (red crosses). 
V2 positions in ye ars 201 0 -2014 a r e show n by empty circles. The black squares indicate the locations of VI in 2012, IBEX RC (assumed 
to represent Bjs , iFunsten et alj (|2009|), downfield), the d irection opposite to RC (upheld), the alte rnative ISMF direction obtained 
from CASSINI/INCA measurements dPialvnas et al .11 2013Tb and the interstellar He llWitte et al .112004 ) and H inflow directions. The He 
inflow direction is taken as the direction of the interstellar flow Vis- The great circle corresponding to the (Bis,Vis) plane (in other 
words He-RC plane) is shown as a red line, the HDP is black (practically coincident with red). The dotted line shows the alternative 
(Bis i Vis) plane corresponding to the assumption that Bis I s given by t he RC position for the highest IBEX energy (triangle) . Blue 
contours show the high pressure region and pressure maximum location IlSchwadron et aU 120141 : iMcComas &; Schwadronl |20 1 4| ). The 
dashed blue line is the great circle (defined by locations of related plasma pressure maximum and the V2 position) on which the crosses 
co rresponding to velocity measure ments should lie if the flow were organized by pressure gradients due to pressure maximum as envisages 
bv lMcComas &; Schwadronl | |20 1 4 ). Thin blue line is the approximate location of the IBEX ribbon. 


with this interpretation. Since V2 position is close to the 
(Bis,Vis) plane, the plasma velocity field measured by V2 
should be approximately parallel to this plane. We find that 
this is indeed the case for the plasma flow vectors which 
come from V2 measurements since 2010.5. 


3 PLASMA VELOCITY OBSERVATIONS 

The V2 trajectory inclination relative to the (Bis , Vis) 
plane has been small from ~ 2005 (Fig. [T] Fig. [5J. At the 
end of 2014 it is still less than 2° (Fig. [2j. Fig. [3] shows 
that, over the last 10 years, excepting abour two years af¬ 
ter the termination shock crossing by V2, the plasma flow 
vectors measured by V2 have been approximately parallel 
to (Bis,Vis) plane, a clear confirmation that this plane is 
indeed an approximate symmetry plane of the flow. 

In Figure |T| the V2 positions and measurements are 
shown projected on the all-sky map. We assume that the 
Bis direct ion is t h e sam e as t he position of the RC ac¬ 
cording to IFunsten et al.l (120091 ) (filled squ are), which is 
close to the energy-averaged RC position llFunsten et al.l 



years 


Figure 2. Inclination angle of V2 radial direction relative to the 
(Bis,Vis) plane versus time. 


12013)1 ) (empty square). The great circles representing the 
(Bis,Vis) plane (red line) and the HDP (black line) are 
then quite close to each other. For measurements between 
2010 and 2014, the V2 positions (empty circles) and the 
yearly averaged plasma velocity vectors (black crosses) are 
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Figure 3. The inclination angle of the proton velocity vectors 
to the (Bis,Vis) plane versus time. Shown are daily and 25- 
day averages (dots and line, respectively). Crosses correspond to 
the angles obtained using t he corrected plasma velocity vectors 
dRichardson Sz Decker|[2014l) . The Bis and V/s are defined as in 
Figure fTI 
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Figure 4. The angle between the proton velocity (measured by 
V2 between 1990 and 2015) and the V2 direction. Shown are daily 
and 25-day averages (dots and line, respectively). Crosses corre¬ 
spond t o the an gles obtained using the corrected plasma velocity 
vectors dRichardson Sz Decked[20141 . 



Figure 5. Blow-up of part of Figure [T] showing the evolution 
of observed plasma velocity direction during V2 journey out¬ 
wards. Th e crosses show the V2 data c orrected for the instrument 
response dRichardson &; Decked 12014) (six-month averages, red) 
and obtained from the raw data (yearly averages, black). The 
data points are numbered (corrected data) or marked by letters 
consecutively in time. The numbers indicate half year periods (1: 
2010.25, 2: 2010.75, 3: 2011.25, 4: 2011.75, 5: 2012.25, 6: 2012.75, 
7: 2013.25, 8: 2013.75, 9: 2014.10). The 9th period is shorter. The 
letters correspond to years (a: 2010, b: 2011. c: 2012, d: 2013, e: 
2014). 


shown. The red crosses show the six-month averages of 
plasma velocity direction which was instrumentally cor¬ 
rected (jRichardson fc Deckerll2014l ). 

It can be seen how well the measured velocity directions 
follow the (Bis Mis) plane being also close to the HDP. For 
comparison, we also show (by the dotted line) the position 
of the (Bis, Vis) plane for another choice of the Bis direc¬ 
tion corresponding to the position_of RC obtai ned from the 
highest energy IBEX data dFunsten et al.|[2013)h (triangle). 
It is clear that the (Bis, Vis) plane defined in this way is 
not a symmetry plane of the flow. 


4 DISCUSSION 

The restriction of the plasma flow vector at V2 to 

the_ (BiSjVis) plane e xplain s also the observation by 

iMcComas fe Schwadronl (120141 ) that, as V2 moves deeper 
into the heliosheath, the flow speed rotates more in the 
transverse than in the polar direction. In Figure[T] this would 
be represented by the flow direction moving along the great 
circle (representing the (Bis,Vis) plane) away from the he¬ 
liospheric nose and from the V2 trajectory direction and 
consequently towards the region where the circle is approx¬ 
imately tangent to the local heliographic parallel. The flow 
direction would therefore have a substantially larger trans¬ 
verse component than the direction of V2 trajectory, while 
staying (approximately) in the same plane as the V2 trajec¬ 
tory and Vis direction. 

The symmetry relative to the (Bis,Vis) plane (de¬ 
termined solely by the parameters of the local interstellar 
medium) must be to some degree broken by the structure of 
the inner heliosphere. In particular, the solar wind is known 
to consist of two components (slow and fast) with the dis¬ 
tribution dependent on the heliographic latitude. However, 
the sola r win d energy flux (ram pressure) is more symmet¬ 
ric dLe Chat et al.|[2012l ), which may reduce the effect on the 
global heliospheric structure. 

The solar magnetic field (in particular, the Parker spi¬ 
ral) is not symmetric relative to the (Bis,Vis) plane. How¬ 
ever, its effect on the bulk plasma flow_in the hchosheath is 
likely to be small. It was shown bv lMcComas fc Schwadron] 
120141) based on separation of the IBEX ENA intensities 
into a part coming from the energetic ions outside the HP 
(IBEX Ribbon) and another part due t o ions in the^in- 
ner h eliosheath (Global Distributed Flux) (ISchwadron et all 
12014) , that the energy density of the heliosheath plasma, de¬ 
termined mainly by the non-thermal ions, exceeds by far the 
magnetic energy density estimated from the VI and V2 data 
(plasma-beta ~20). This means that the magnetic field in 
the heliosheath can not induce a significant departure from 
the sym metr y imp ressed b y the external situation. 

IMcComas fe Schwadronl 120141 ) relate the V2 measure¬ 
ments of the heliospheric plasma flow to distribution of the 
integrated plas ma pressure infer red from IBEX ENA obser¬ 
vations IlSchwadron et all 120141 ). They identify the region 
of high lieliosheath plasma pressure (approximately repre¬ 
sented by a blue contour in our Figure [lj and argue that 
the plasma flow is directed away from the point of max¬ 
imum pressure (blue circle in our Figure [TJ, which differs 
from the position of the heliospheric nose (He inflow di¬ 
rection). However, we note that, in this scenario, the flow 
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velocity direction would lie in the plane defined by the V2 
position and the point of maximum pressure (dashed blue 
line in our Figure [lj. 

As shown in Figure [T| V2 measurements follow the 
(Bis,Vis) plane. Combination of this (approximate) sym¬ 
metry with the geometrical coincidence of V2 position (in 
years 2013-2015) with the (Bis,Vis) plane explains there¬ 
fore in a natural way the averaged observed plasma ve¬ 
locity vectors (crosses in Figure [TJ without recurring to 
the idea of heliosheath plasma being pushed away from 
th e excess pressure regi on near the nose as speculated in 
iMcComas fc Schwadronjjl gfn^db_ 

Nevertheless, McComas fe Schwadronl (12014) idea that 
the flow should be linked with the pressure distribution 
is important. If a strict mirror symmetry relative to the 
(Bis,Vis) plane were valid, one would expect that the high 
plasma pressure region should have the same symmetry. Fig¬ 
ure [l] shows that the high pressure region is bisected by the 
(Bis,Vis) plane, but not into equal parts. However, the 
edge of the high pressure region near the V2 position is 
approximately perpendicular to the (Bis,Vis) plane. This 
implies that the local pressure gradient is likely to cause the 
plasma flow along the plane, as observed by V2. 

The asy mmetry o f _the pres sure distribution 
obt ained by Schwadron et_al .1 (120141 ^1 and used by 

IMcComas fc Schwadronl ( 20141 ') is not necessarily in 
conflict with our results. This is because the pressure 
they consider is integrated over the line-of-sight (LOS). 
It is therefore affected also by the region close to the 
termination shock, where V2 velocity measurements deviate 
significantly from the (Bis, Vis) plane (see Figures 03 0. 
The velocity measured during first half of the year 2010 is in 
fact close to the dashed blue line in Figures [T] and |3j which, 
as we have pointed o ut above, corresponds to the pr essure 
distribution used by IMcComas fe Schwadronl (1201411 . The 
symmetry with respect to the (Bis, Vis) may then hold 
for the local (not integrated) pressure distribution at larger 
distances from the termination shock, where V2 velocity 
measurements follow the (Bis, Vis) plane. 


5 CONCLUSIONS 

V2 position near the HDP plane permits a check of mir¬ 
ror symmetry of the heliosheath (inner and outer) plasma 
flow with respect to this plane. Our main result is that V2 
observations are consistent with this symmetry. 

Assuming that the Bis direction is given by the IBEX 
RC and the Vis direction by the the neutral interstellar 
He flow, the plasma velocity measured by V2 lies in the 
(Bis,Vis) plane (almost identical with the HDP). As V2 
approaches the heliopause, the direction of plasma velocity 
rotates, still within the (Bis,Vis) plane. This explains why 
the transverse component of the plasma velocity at V2 be¬ 
comes larger than the polar one. 

The V2 data suggest that the mirror symmetry with 
respect to (Bis,Vis) plane is restricted to the outer part of 
the inner heliosheath. The plasma pressure distribution ob- 
tained from IBEX data is integrated over the full LOS length 
l|Schwadron et ahll2014l h This may explain why it departs 
from this symmetry. 

Because V2 location is approximately in the symmetry 


plane, we expect that the heliospheric plasma flow at the 
heliopause will be approximately parallel to the interstellar 
plasma flow just outside the heliopause, and anti-parallel to 
the interstellar magnetic held. Figures [T] [4] and [5] show that 
the plasma velocity at V2 location in fact evolves towards 
the anti-held direction. 

At the heliopause, the draped ISMF would, in general, 
differ from the undisturbed ISMF. However, V2 trajectory 
directio n is close to perpen dic ular ( 97° ) relat iv e to the 
Bis ■ In ICrygorczuk, Czechowski fe Grzedzielski! (l2014f) we 
showed that the magnetic held immediately outside the he- 
liopause observed by V2 should then be similar in direction 
to the undisturbed ISMF (RC). 
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